(19) 



J 



Europaisches Patentamt 
European Patent Office 
Office europeen des brevets 



(12) 



(11) EP 1 275 967 A1 

EUROPEAN PATENT APPLICATION 



(43) Date of publication: 

15.01.2003 Bulletin 2003/03 

(21 ) Application number: 02077641 .5 

(22) Date of filing: 17.12.1998 



(51) intCJT: G01 P 5/1 8, G01 F 1/708, 
G01P 5/10 



(84) Designated Contracting States: 


• Schuldt, Spencer B. 


DE ES Fl FR GB IT NL SE 


Bloomingtom, Minnesota 55437 (US) 




• Ulrich, Bonne 


(30) Priority: 31.12.1997 US 2157 


Hopkins, Minnesota 55343 (US) 


(62) Document number(s) of the earlier application(s) in 


(74) Representative: 


accordance with Art. 76 EPC: 


Fox-Male, Nicholas Vincent Humbert et al 


96964002.4 / 1 049 935 


Eric Potter Clarkson 




Park View House 


(71) Applicant: Honeywell Inc. 


58 The Ropewalk 


Minneapolis Minnesota 55440-0524 (US) 


Nottingham NG1 5DD (GB) 


(72) Inventors: 


Remarks: 


• Kubisiak, David 


This application was filed on 03 - 07 - 2002 as a 


Chanhassen, Minnesota 55317 (US) 


divisional application to the application mentioned 


• Matthys, Robert J. 


under INID code 62. 


St. Anthony, Minnesota 55418 (US) 





(54) Device and method for measuring fluid velocity by determinig the time-of-flight of heat 
pulses 



(57) A method and apparatus for measuring the ve- 
locity of a fluid relatively independently of the physical 
properties of the fluid. This is preferably accomplished 
by spacing two sensor elements at different distances 
from the heater element. The present invention also 
contemplates minimizing the effects of the non-zero 
heater time lag and/or the non-zero sensor time lag for 



increased accuracy. This is preferably accomplished by 
either measuring the time lags and subtracting the val- 
ues thereof from an uncorrected transit time measure- 
ment, forcing the sensor elements to track the thermal 
disturbance in the fluid thereby minimizing the effects of 
the sensortime lags, or measuring the transit time using 
sensors that have substantially zero thermal mass. 
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Description 

BACKGROUND OF THE INVENTION 

5 1. FIELD OF THE INVENTION 

[0001] The present invention generally relates to the measurement of fluid properties and, more particularly, to the 
determination of the velocity of a fluid of interest. 

10 2. DESCRIPTION OF THE PRIOR ART 

[0002] A number of techniques have been devised to measure the velocity of a fluid within a lumen. One approach 
is the "time of flight" approach, which generally includes determining the time require for a thermal wave to flow from 
a source heater element to a destination sensor element. By knowing the distance between the heater and sensor, 

15 the velocity of the fluid can thus be calculated. 

[0003] U.S. Patent No. 4,576,050 to Lambert discloses one such 'lime of flight approach. Lambert suggests ener- 
gizing a heater strip with an oscillating heater input signal to emit thermal waves in the fluid. The thermal waves prop- 
agate through the fluid at a rate that is dependent on the fluid velocity that flows perpendicular to the heater strip. A 
thermo-electric detector, spaced from one or both side of the heater, senses the thermal wave and provides a corre- 

20 sponding detector output signal. The velocity of the fluid can be determined, at least to first order, from the time differ- 
ential between the heater input signal and the detector output signal. 

[0004] A limitation of the Lambert approach is that the measured time or phase differential between the heater input 
signal and the detector output signal depends on a number of physical properties of the fluid, including, for example, 
temperature, pressure, thermal conductivity, and thermal diffusivity. To compensate for some of these parameters, 

25 Lambert suggests making a reference measurement with a second sensor that is exposed to still fluid or fluid having 
a velocity component that is perpendicular to the source of the second sensor. By comparing the output of the second 
sensor with that of the first sensor, Lambert suggests that a phase difference can be computed that is independent of 
at least some of the fluid properties listed above. As can be seen, the second sensor of Lambert is not used to measure 
the fluid velocity, but rather is used to compensate for selected physical properties of the fluid. 

30 [0005] Another limitation of the Lambert approach is that several potential error sources in the phase differential 
measurement are neglected, thereby reducing the accuracy of the measurement. One such error source is the non- 
zero heater time lag that typically exists between the heater input signal and the elevated temperature response of the 
heater element (and thus the fluid). The heater time lag is typically dominated by the thermal conductivity, k, of the 
fluid of interest, at least for microbridge structures as contemplated by a preferred embodiment of the present invention. 

35 [0006] Another error source is the non-zero sensor time lag that typically exists between the arrival of the temperature 
disturbance at the sensor element and the corresponding response of the sensor element. The temperature of the 
sensor element typically does not react instantaneously to a temperature change in the fluid, primarily due to the 
sensors non-zero thermal mass. 

[0007] It would be desirable, therefore, to provide an method and apparatus whereby the fluid velocity can be de- 
40 termined relatively independently of the physical properties of the fluid. It would also be desirable to provide a flow 
sensor that accounts for the non-zero heater time lag and the non-zero sensor time lag for increased accuracy. 

SUMMARY OF THE INVENTION 

45 [0008] The present invention provides apparatus as defined in Claim 1 hereinafter. 
[0009] The apparatus may include the features of dependent Claim 2. 
[0010] The present invention also provides a method as defined in Claim 3. 
[0011] The method may include the features of Claim 4. 

[0012] The present invention may overcome many of the disadvantages associated with the prior art by providing a 
so method and apparatus whereby the fluid velocity can be determined relatively independently of the physical properties 
of the fluid, and minimizes the effects of the non-zero heater time lag and/or the non-zero sensor time lag for increased 
accuracy. 

[001 3] In a first embodiment of the present invention, a heater element and at least two spaced sensor elements are 
provided, whereby the sensor elements are spaced at different distances from the heater element. The heater and 
55 sensor elements are disposed in and closely coupled to a fluid medium (fluid or gas) of interest. A time-varying input 
signal is provided to the heater element, which heats the surrounding fluid. Of interest is the transit times for the 
temperature disturbance to travel from the heater element to selected sensor elements. Because the sensor elements 
are spaced at different distances from the heater element, and as more fully described below, the fluid velocity can be 
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calculated from the measured transit time values relatively independently of the fluid properties. 
[0014] To reduce the effects of the non-zero heater and sensor time lags on the desired transit times, a first embod- 
iment of the present invention contemplates measuring the heater and sensor time lags during a calibration procedure. 
The heater and sensor time lags can be directly measured by providing a heater input signal to the heater element, 

5 and a sensor input signal to each of the sensor elements. By monitoring the transient temperature responses of the 
heater and sensor elements via heater and sensor output signals, respectively, the heater and sensor time lags can 
be determined. To obtain more accurate transit times for the temperature disturbance in the fluid, the heater time lag 
and/or the sensor time lags may be subtracted from the overall time lags measured between the heater input signal 
and each of the sensor output signals. 

w [0015] Another approach to reduce the effects of the non-zero heater and sensor time lags is to provide heat to the 
sensor elements via the corresponding sensor input signal during operation. The sensor input signals are preferably 
controlled to provide a frequency, phase and amplitude that produce a resulting temperature response in the corre- 
sponding sensor element that tracks the temperature disturbance in the fluid. When this condition is satisfied, substan- 
tially zero heat is transferred from the fluid to the sensor elements, and the effect of the sensor time lags are minimized. 

is [0016] The proper phase and amplitude of the sensor input signals can be determined during a calibration procedure. 
In one embodiment, the sensor elements are first subjected to a vacuum condition, and a sensor input signal is provided 
to each of the sensor elements. Because no fluid surrounds the sensor elements, substantially no heat is transferred 
from the sensor elements to the fluid. A sensor time lag between the sensor input signal and the resulting temperature 
response of each sensor element is measured and stored. A power/resistance ratio of each sensor element is also 

20 measured and stored. 

[0017] The sensor elements are then subjected to a fluid of interest. During operation, the phase of each of the 
sensor input signals is adjusted so that the resulting sensor time iag equals the sensor time lag measured under vacuum 
conditions. Likewise, the amplitude of each of the sensor input signals is adjusted so that the resulting power/resistance 
ratio equals the powcr/rcsistancc ratio measured under vacuum conditions. Under these conditions, the sensor time 

*5 lags may be reduced or eliminated as a factor when determining the transit times of the temperature disturbance. - 
[001 8] The effect of the heater lag time may be reduced or eliminated by beginning the time lag measurement when 
a heater output signal crosses a predetermined threshold, where the heater output signal is proportional to the resist- 
ance (temperature) of the heater element. Alternatively, the effect of the heater time lag may be reduced or eliminated ; 
by starting the time lag measurement when the heater input signal crosses a predetermined threshold, and then sub=- 

30 trading the heater time lag from the resulting overall time lag value. 

[0019] Rather than measuring the temperature disturbance with a temperature sensitive resistor, it is contemplated 
that the velocity of the fluid may be determined by using optical sensors. In an illustrative embodiment, a time lag is • 
determined between a first location and a second location for a transient elevated temperature condition in the fluid., 
using a first and second optical sensors. The first optical sensor optically detects the transient elevated temperature 

35 condition at the first location . and the second optical sensor optical detects the transient elevated temperature condition; 
at the second location. By knowing the distance between the first and second optical sensors, the fluid velocity may:* 
be determined. It is also contemplated that optical sensors may be used to measure the various lag times to determine 
the thermal conductivity, thermal diffusivity and specific heat of a fluid of interest. 

40 BRIEF DESCRIPTION OF THE DRAWINGS 

[0020] Other objects of the present invention and many of the attendant advantages of the present invention will be 
readily appreciated as the same becomes better understood by reference to the following detailed description when 
considered in connection with the accompanying drawings, in which like reference numerals designate like parts 
45 throughout the figures thereof and wherein: 



Figures 1, 2 and 3 are different views of a prior art embodiment of a microbridge flow sensor; 
Figure 4 is a partial cut-away view of a microbridge sensor system; 

Figure 5 is a first illustrative cross sectional view of a microbridge sensor in accordance with the present invention 
50 having two downstream sensor elements; 

Figure 6 is a second illustrative cross sectional view of a microbridge sensor in accordance with the present 
invention having an upstream and a downstream sensor clement; 

Figure 7 is a third illustrative cross sectional view of a microbridge sensor in accordance with the present invention 
having more than two downstream sensor elements; 

Figure 8 is a fourth illustrative cross sectional view of a microbridge sensor in accordance with the present invention 
having a number of downstream and upstream sensor elements; 

Figure 9 is a timing diagram showing the desired transit time Az t , the heater lag time Az h , and the sensor lag time 

AZ.; 



55 



DCCID: <EP . 1275967A1_I..> 



3 



EP 1 275 967 A1 



Figure 10 is a schematic diagram of an illustrative circuit for use with the microbridge heater and sensor elements 
of, for example, Figure 5; 

Figure 11 A and Figure 11B are a block diagram showing another illustrative embodiment of the present invention 
for obtaining the various time lags that are used to determine the velocity, the thermal conductivity, thermal diffu- 
5 sivity, and specific heat of the fluid of interest; 

Figure 12 is a timing diagram showing an illustrative input signal, and an inverted copy thereof, which are provided 
by the signal generators shown in Figure 11A and Figure 11 B; 

Figure 13 is a timing diagram showing an illustrative high frequency signal modulated by the input signals of Figure 
12; 

10 Figure 14 is a timing diagram showing the power of the modulated input signal of Figure 13; 

Figure 15 is a timing diagram showing the resistance of the heater and/or sensor elements versus time when the 
power signal of Figure 14 is provided thereto; 

Figure 16 is a timing diagram showing the sensor time lag Az^^ when measured under vacuum conditions; 
Figure 17 is a timing diagram showing the various time lags measured by the embodiment of Figure 11A and 
15 Figure 11 B; 

Figure 18 is a block diagram showing a processing means that is used for calculating the velocity, thermal con- 
ductivity and specific heat of the fluid of interest using the time lag values provided by Figure 1 1 A and Figure 1 1 B; 
Figure 19 is a schematic diagram of an illustrative embodiment that uses optical sensor means to measure the 
transit time of a thermal disturbance from a first location to a second location; and 
20 Figure 20 is a schematic diagram of another illustrative embodiment that uses optical sensor means to measure 

the transit time of a thermal disturbance from a first location to a second location. 



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

25 [0021] The present invention, then, is directed to a system that determines the fluid velocity relatively independently 
of the physical properties of the fluid, and minimizes the effects of the non-zero heater time lag and/or the non-zero 
sensor time lag for increased accuracy. The preferred embodiments of the present invention contemplate disposing a 
microscopic sized heating element in a sample of the fluid of interest. The microsensor system or "microbridge", as it 
will be referred to herein, though not limiting, is presently preferred for several reasons. The system is extremely fast- 

30 reacting, is very accurate, very sensitive because of its advantageous coupling to the fluid of interest and small and 
adaptable to a variety of configurations. 

[0022] The microbridge semiconductor chip sensor contemplated, for example, in certain embodiments preferred 
for the invention may resemble the form of one or more of the microbridge systems illustrated in U.S. Patent No. 
4,478,076, U.S. Patent No. 4,478,077, U.S. Patent No. 4,501,144, U.S. Patent No. 4,651,564, U.S. Patent No. 

35 4,683,159, and U.S. Patent No. 4,994,035, all of common assignee with the present invention. 

[0023] Such a system is exemplified by Figures 1-3 taken from U.S. Patent No. 4,994,035 to Aagard et al. A dis- 
cussion of that example will now be presented as it will be helpful in understanding the present invention. While the 
present discussion is believed sufficient, to the extent necessary, any additional material contained in the microbridge 
related patents cited is deemed to be incorporated herein by reference. 

40 [0024] The prior art system of Figures 1-3 contemplates a pair of thin film temperature sensors 22 and 24, a thin 
film heater 26 and a support member 20 supporting the sensors and heater out of contact with the base substrate. 
Sensors 22 and 24 are disposed on opposite sides of heater 26. Support member 20 is a semiconductor, preferably 
silicon, chosen because of its adaptability to precision etching techniques and ease of electronic chip producibility. The 
embodiment includes two identical temperature sensing resistor grids 22 and 24 acting as the thin film heat sensors 

45 and a centrally located heater resistor grid 26 acting as the thin film heater. 

[0025] Sensors 22 and 24 and heater 26 may be fabricated of any suitable, stable metal or alloy film. The metal used 
may be a nickel-iron alloy sometimes referred lo as permalloy, wilh a composition of 80 percent nickel and 20 percent 
iron. The sensor and heater grids are encapsulated in a thin film of dielectric, typically comprising layers 28 and 29 
and preferably silicon nitride, Si 3 N 4 to form the film members. 

50 [0026] In Figures 1 and 2, the sensor comprises two thin film members 32 and 34, with member 32 comprising 
sensor 22 and member 34 comprising sensor 24, each member comprising one-half of heater 26 and having a preferred 
dimension of 150 microns wide and 400 microns long. 

[0027] The system further describes an accurately defined fluid space 30 that effectively surrounds elements 22, 24, 
26, and is achieved by fabricating the structure on silicon surface 36. Thin film elements 22, 24 and 26 have thicknesses 
55 of approximately 0.08 to 0.12 micron with line widths on the order to 5 microns and spaces between lines on the order 
of 5 microns. The elements encapsulated in the silicon nitride film preferably have a total thickness of approximately 
0.8 microns or less. The fluid space 30 may be fabricated by subsequently etching an accurately defined fluid space 
of about 1 00 microns deep into silicon body 20 beneath members 32 and 34. 
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25 



[0028] Members 32 and 34 connect to top surface 36 of semiconductor body 20 at one or more edges of depression 
or fluid space 30. As illustrated in Figure 3, members 32 and 34 may be bridged across depression 30; alternately, for 
example, members 32 and 34 could be cantilevered over depression 30. 

[0029] In the system shown, heat flows from the heater to the sensor by means of both solid and fluid couplings 
therebetween . Of note is the fact that silicon nitride (Si 3 N 4 ), besides being a good electrical insulator, is also an effective 
solid thermal insulator. Because the connecting silicon nitride film within members 32 and 34 is a good insulator, heat 
transmission through the solid does not dominate the propagation of heat from heater 26. This further enhances the 
relative amount of the heat conducted to sensing resistors 22 and 24 from heater resistor 26 by flow through the 
surrounding fluid rather than through the supporting nitride film. Moreover, the supporting silicon nitride film has a low 
enough thermal conductivity that sensing resistor grids 22 and 24 can be located immediately adjacent or juxtaposed 
to heating resistor grid 26. Thus, sensing resistor grids 22 and 24 are in effect suspended rigidly in the fluid space 
proximate heater resistor 26 and act as thermal probes to measure the temperature of the air near and in the plane of 
heater resistor grid 26. 

[0030] Figure 4 is partial cut-away view of a microbridge sensor system placed in line with a flow pipe. A main flow 
'5 channel 200 having a central bore 202 is connected to the pipe that carries a fluid of interest. A first chamber 204 is 
in fluid communication with the central bore 202 of the flange 200 via a single lumen 206. A header 208 having a first 
microbridge sensor 21 0 mounted thereto is inserted into the first chamber 204 and secured to the main flow channel 
200 as shown. In this configuration, the first microbridge sensor is exposed to the fluid of interest with substantially 
zero flow. The first microbridge sensor 210 is typically used to measure fluid properties such as thermal conductivity, 
20 thermal diffusivity, specific heat, temperature and pressure. 

[0031 ] A second sensor 222 is positioned in a small bypass channel 21 4. In this configuration , the second microbridge 
sensor is exposed to the flow of the fluid of interest. The second microbridge sensor 222 is typically used to measure 
fluid velocity. 

[0032] Figure 5 is a first illustrative cross sectional view of a microbridge sensor in accordance with the present 
invention, having two downstream sensor elements. As indicated above, the present invention provides apparatus and 
methods for determining the fluid velocity relatively independently of the physical properties of the fluid, and may also 
minimize the effects of the non-zero heater time lag and/or the non-zero sensor time lag for increased accuracy. 
[0033] Generally, the present invention includes a heater element 224 and at least two spaced sensor elements .226 
and 228. The sensor elements 226 and 228 are preferably spaced different distances from the heater element 224. In 
30 the illustrative diagram, sensor 226 is spaced a first distance "ay* from the heater element 224, and sensor 228 is 
spaced a second distance "d^* from the heater element 224. Both sensor 226 and 228 are shown downstream from 
the heater element 224. 

[0034] The heater element 224, sensor element 226 and sensor element 228 are preferably disposed in and closely- 
coupled to a fluid medium (fluid or gas) of interest. A time-varying input signal is provided to the heater element 224,. 
35 which heats the surrounding fluid. Of interest is the transit time for the temperature disturbance in the fluid to travel . 
from the heater element 224 to each of the sensor elements 226 and 228. The fluid velocity can be calculated from- 
this transit time value, relatively independently from the fluid properties. 

[0035] The heater element 224 is shown having a support member 228 that supports the heater element 230 out of 
contact with the base substrate 232. Together, the heater element 224 and support member 230 form a heater film 

40 member. Likewise, the sensor element 226 is shown having a support member 234 that supports the sensor element 
226 out of contact with the base substrate 230. Together, the sensor element 226 and support member 234 form a first 
sensor film member. Finally, the sensor element 228 is shown having a support member 236 that supports the sensor 
element 228 out of contact with the base substrate 230. Together, the sensor element 228 and support member 236 
form a second sensor film member. 

45 [0036] Heater element 224 and sensor elements 226 and 228 may be fabricated of any suitable, stable metal or 
alloy such as platinum, Nickel, Iron-Nickel, etc. Heater element 224 and sensor elements 226 and 228 may be any 
resistive element including a wire, but are preferably a film. Moreover, healer element 224 and sensor elements 226 
and 228 may be of any shape including a grid pattern as described above, or simply a line. As indicated above, the 
heater element 224 and sensor elements 226 and 228 are preferably encapsulated in a thin film of dielectric, such as 

so silicon nitride, Si 3 N 4 , to form the support members 230, 234 and 236. 

[0037] An accurately defined fluid (gas or liquid) space 240 is preferably provided which effectively surrounds heater 
clement 224 and sensor elements 226 and 228, and is achieved by fabricating the structure on silicon surface 242. 
Heater element 224 and sensor elements 226 and 228 preferably have thicknesses of approximately 0.08 to 0.12 
micron, with line widths on the order to 5 microns and, if a grid is used, spaces between lines on the order of 5 microns. 

55 The fluid space 240 may be fabricated by subsequently etching an accurately defined fluid space of about 1 00 microns 
deep into silicon substrate 232 beneath heater element 224 and sensor elements 226 and 228. 
[0038] Support member 230 and heater element 224 preferably connect to top surface 242 of semiconductor sub- 
strate 232 at one or more edges of depression or fluid space 240. Support member 230 and heater element 224 may 
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bridge across depression 240 as shown, or alternately, for example, cantilevered over depression 240. The sensor 
elements 234 226 and 228 are preferably similarly constructed. It is recognized that any number of heater and sensor 
elements may be provided in a like manner. However, for illustration purposes, only one heater element 224 and two 
sensor elements 226 and 228 are shown in Figure 5. 

5 [0039] A heater input signal is preferably provided to the heater element 224, which induces a thermal disturbance 
in the fluid. Each of the sensor elements 226 and 228 sense the arrival of the thermal disturbance at their respective 
locations. Of interest is the transit times for the temperature disturbance to travel from the heater element 224 to each 
of the sensor elements 226 and 228. As described more fully below, because the sensor elements 226 and 228 are 
spaced at different distances from the heater element, the fluid velocity can be calculated relatively independently of 

10 the fluid properties. 

[0040] Before providing the relation for fluid velocity as a function of the transit time values between the heater 
element and two sensor elements, namely Az 1 and Az 2 , some background information is provided. The propagation 
or spread of a local, non-uniform temperature distribution can be derived from the following three-dimensional temper- 
ature diffusivity or heat conductivity equation, to which a term has been added to represent fluid velocity in the x- 
*s direction: 

dT/dz = D t 9 2 T/ax 2 - V x 8T/dx (1) 
20 The one dimensional transient solution to equation (1 ) is: 

T ~ z^expHdj - vz) 2 /(4D t z)) (2) 

25 where V has been replaced by (± dj ± vz) to include the fluid velocity term, and dj represents the distance from the 
heater to the "i thM sensor. To determine the time lag between very short forcing functions at time z=0 and the arrival of 
these pulses in terms of maxima in T(z) at the sensor element position dj, equation (2) is differentiated and set to zero: 

30 dT/dz = 0 - -1/z 2 + (1/z)(2v(dj - vz)/4D t z + (dj - vz) 2 /4D t z 2 ) (3) 

multiplying by 7 2 and rearranging yields the relation: 

35 -4D t z - v 2 z 2 + dj 2 = 0 (4) 

As shown by equation (4), choosing d 1 = -d 2 does not cause a change in the time lag because there is no linear term in d. 
[0041] To determine a relation that is independent of the fluid properties, which are here represented by D t in equation 
^ (4), two unequal values of dj are considered, namely d 1 and d 2 (d 1 1 d 2 ). This yields: 

v% = <i*/z A -4D t (5) 
45 v 2 z 2 = d 2 2 /z 2 - 4D t (6) 

v 2 = (d*/z, - d 2 2 /z 2 )/( Zl - z 2 ) (7) 



50 



55 



Equation (7) has been found to accurately represent the fluid velocity over a wide range of fluid velocities from, for 
example, 1 6,384 cm/s down to 32 cm/s, relatively independently of the physical properties of the fluid. It is recognized 
that equation (7) may not model the physics at very low flow rates. At low flow rates, selected physical properties such 
as thermal diffusivity may become a factor. It is contemplated, however, that by properly locating the sensor elements, 
these effects can be minimized even at low flow rates, as further described below. 

[0042] Rather than providing both sensors downstream from the heater element as shown in Figure 5, it is contem- 
plated that one sensor element 250 may be placed upstream and another sensor 252 may be placed downstream of 
the heater element 254, as shown in Figure 6. To use equation (7), sensor 250 is spaced a different distance from the 
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heater element 254 than sensor 252. 

[0043] To reduce the possible negative effects of thermal diffusivity and other properties of the fluid at low flow rates, 
it is contemplated that a first set of sensor elements may be used for measuring low flow rates and another set may 
be used for higher flow rates. For example, in Figure 7, those sensor that are positioned closest to the heater element, 
5 such as sensor elements 280 and 282, may be used to measure higher flow rates, as the thermal diffusivity component 
may be negligible at the higher flow rates. Likewise, sensor elements that are positioned further from the heater element 
may be used to measure the lower flow rates, including sensor 284. Using this approach, the effect of the thermal 
diffusivity component on the flow rate measurement may be minimized. 

[0044] In addition, it is contemplated that a higher amplitude heater input signal may be provided when measuring 
io high flow rates, and conversely, a lower amplitude heater input signal may be provided when measuring low flow rates. 
A higher amplitude temperature disturbance can be more easily detected, but can increases the speed of the thermal 
diffusivity component in the fluid. Thus, a lower amplitude heater input signal may reduce the speed of the thermal 
diffusivity component, and provide more accurate results at lower flow rates. 

[0045] Figure 8 is a fourth illustrative cross sectional view of a microbridge sensor in accordance with the present 

is invention having a number of downstream and upstream sensor elements. In this embodiment, a number of sensor 
elements are equally spaced from the heater element both in an upstream and downstream direction. In operation, 
however, only the outputs of selected sensor elements that are spaced at different distances from the heater element 
are selected, and thus the fluid velocity may be obtained using equation (7). For example, the closest sensor element 
300 in the downstream direction may be used in conjunction with the second closest sensor element 302 in the upstream 

20 direction. Other combinations are also contemplated. 

[0046] One approach to measuring the transit time of the temperature disturbance from the heater element to a 
sensor element is to start a timer when the heater, input signal crosses a predetermined threshold, and stop the timer 
when the temperature response of the sensor element crosses a predetermined threshold. As described above, the 
fluid velocity can bo determined by using the transit time values for two or more sensors that arc spaced at different 

25 distances from the heater element, relatively independent from the fluid properties. 

[0047], While such an approach may be sufficient for many applications, several factors should be considered for 
higher precision applications. A first factor is the non-zero heater time lag Az h between the heater input signal and the 
elevated temperature response of the heater element (and fluid). An illustrative heater input signal and elevated tem- 
perature response of the heater element are shown at 320 and 322 of Figure 9, respectively. The heater time lag Az h 

30 is typically dominated by the thermal conductivity, k, of the fluid of interest, at least for microbridge structures as con- 
templated by a preferred embodiment of the present invention. As discussed in U.S. Patent Application Serial Number. 
09/002,156, entitled "METHOD AND APPARATUS FOR MEASURING SELECTED PROPERTIES OF A FLUID OF 
INTEREST USING A SINGLE HEATER ELEMENT*, the thermal conductivity, k, of the fluid of interest can be derived 
from the heater time lag, when measured with substantially zero fluid flow. 

35 [0048] A second factor is the non-zero sensor time lag Az s between the arrival of the temperature disturbance at the 
sensor element and the resistance change in the sensor element. That is, the temperature of the sensor element 
typically does not react instantaneously to a temperature change in the fluid, primarily because of the sensors non- 
zero thermal mass. This is explicitly shown in Figure 9 wherein the arrival of the temperature disturbance at the sensor 
is shown at 324, and the resistive response of the sensor element is shown at 326. 

40 [0049] One approach to reduce the negative impact of Az h and Az s is to measure the heater and sensor time lags 
during a calibration procedure. This can be accomplished by providing a heater input signal to the heater element and 
a sensor input signal to each of the sensor elements under flow conditions. The heater time lag Az h and each of the 
sensor time lags Az s can then be determined by monitoring the transient temperature responses of the heater and 
sensor elements, respectively. To obtain a more accurate transit time Az t for the temperature disturbance in the fluid, 

45 the heater time lag Az h and the sensor time lag Az s may then be subtracted from the time lag measured between the 
heater input signal and the temperature response of the corresponding sensor element. 

[0050] Another approach is to provide heat to the sensor elements via sensor input signals during operation. The 
sensor input signals are preferably controlled to provide a frequency, phase and amplitude that produce a resulting 
temperature response in the corresponding sensor elements that tracks the temperature disturbance of the fluid. When 
5 ° this condition is satisfied, substantially zero heat is transferred from the fluid to the sensor element. 

[0051] Figure 10 is a schematic diagram of an illustrative circuit for use with the microbridge heater and sensor 
elements of, for example, Figure 5. In accordance with the present invention, a periodic time-varying heater input 
signal 350 is provided to a heater element 224. The heater element 224 is disposed in and closely coupled to a fluid 
medium (gas or liquid) of interest. To simultaneously provide power to, and measure the resistance response of the 
heater element 224, the heater element 224 may be incorporated into one leg of a Wheatstone bridge 352, as shown. 
The time varying heater input signal 350 is provided to a power input terminal 354 of the Wheatstone bridge 352, which 
is either directly or indirectly coupled to the heater element 224. In this configuration, the Wheatstone bridge 352 
provides a differential output signal that has an amplitude that is proportional to the resistance of the heater element 
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224. Preferably, the differential output signal is provided to a differential amplifier circuit 356 to provide an amplified 
output signal 358. 

[0052] The heater input signal 350 provides power to the heater element 224 and induces a transient elevated tem- 
perature condition in the heater element 224 and the fluid of interest. Because the heater element 224 is closely coupled 
to the fluid, the thermal conductivity "k" of the fluid directly effect the time variable temperature response of the heater 
element 224. Further, the thermal conductivity of the fluid is typically dependent on the pressure and/or temperature 
of the fluid. Thus, it has been found that the thermal conductivity, pressure and/or temperature of the fluid of interest 
can be determined by examining a variable phase lag or time lag between the input signal 350 provided to the heater 
element 224 and a subsequent transient temperature response of the heater element 224 when measured with sub- 
stantially zero fluid flow. 

[0053] Likewise, a periodic time-varying sensor input signal 374 may be provided to a first sensor element 226. The 
first sensor element 226 is spaced from the heater element 224 and disposed in and closely coupled to a fluid medium 
(gas or liquid) of interest. To simultaneously provide power to, and measure the resistance response of the first sensor 
element 226, the first sensor element 226 may be incorporated into one leg of a Wheatstone bridge 378, as discussed 
above. The sensor input signal 374 is provided to a power input terminal 380 of the Wheatstone bridge 378, which is 
either directly or indirectly coupled to the first sensor element 226 as shown. In this configuration, the Wheatstone 
bridge 378 provides a differential output signal that has an amplitude that is proportional to the resistance of the first 
sensor element 226. Preferably, the differential output signal is provided to a differential amplifier circuit 382 to provide 
an amplified output signal 384. 

[0054] Finally, it is contemplated that a periodic time-varying sensor input signal 400 may be provided to a second 
sensor element 228. The second sensor element 228 is spaced a different distance from the heater element 224 than 
the first sensor element 226. To simultaneously provide power to, and measure the resistance response of the second 
sensor element 228, the second sensor element 228 may be incorporated into one leg of a Wheatstone bridge 402, 
as discussed above. The sensor input signal 400 is provided to a power input terminal 401 of the Wheatstone bridge 
402, which is either directly or indirectly coupled to the second sensor element 228 as shown. In this configuration, the 
Wheatstone bridge 402 provides a differential output signal that has an amplitude that is proportional to the resistance 
of the second sensor element 228. Preferably, the differential output signal is provided to a differential amplifier circuit 
404 to provide an amplified output signal 406. 

[0055] Using the circuit shown in Figure 10, the heater time lag Az h and the sensor time lags A7 512 can be determined 
by monitoring the transient temperature responses of the heater and sensor elements, respectively. To obtain a more 
accurate transit time Az t for the temperature disturbance in the fluid during operation, the heater time lag Az h and the 
corresponding sensor time lags Az s1 2 may be subtracted from the corresponding time lags measured between the 
heater input signal and the temperature response of the respective sensor elements. 

[0056] Moreover, the circuit shown in Figure 10 may be used to introduce heat into the sensor elements via the 
sensor input signals 374 and 400. The sensor input signals 374 and 400 are preferably controlled to provide a frequency, 
phase and amplitude that produce a resulting temperature response in the sensor elements 226 and 228 that tracks 
the temperature disturbance of the fluid. When this condition is satisfied, substantially zero heat is transferred from the 
fluid to the sensor elements 226 and 228. 

[0057] Figure 1 1 is a block diagram showing another illustrative embodiment of the present invention for obtaining 
various time lags values to determine the velocity, thermal conductivity, thermal diffusivity, and specific heat of a fluid 
of interest. In this embodiment, heat is introduced into the sensor elements via a corresponding sensor input signal 
during operation. The sensor input signals are controlled to provide a frequency, phase and amplitude that produce a 
resulting temperature response in the sensor elements that tracks the temperature disturbance of the fluid. 
[0058] In this illustrative embodiment, both the heater input signal and the sensor input signals include a high fre- 
quency component and a lower frequency component. The lower frequency component preferably modulates the high 
frequency component. This construction allows an accurate dosage of power to be delivered to the heater element 
and Lhe sensor elements, while eliminating frequency doubling effects that may occur without adding a DC component 
to the input signal. The high frequency component is preferably in the range of 0.1 to 3 MHZ, while the lower frequency 
component is preferably in the range of 30 to 200 Hz. 

[0059] With respect to the heater element, a first generator 500 generates and provides the high frequency compo- 
nent of the input signal to a modulator 502 via interface 504. The generator 500 also generates the lower frequency 
component, and an inverted copy thereof, and provides those signals to modulator 502 via interfaces 506 and 508, 
respectively. An illustrative lower frequency component 540 and inverted copy thereof 542 are shown in Figure 12. 
The modulator 502 modulates the high frequency component using the lower frequency component signals to produce 
a modulated heater input signal. An illustrative modulated heater input signal is shown in Figure 13. The power deliv- 
ered by the modulated heater input signal is shown in Figure 14. 

[0060] The modulated heater input signal is provided to heater block 510. To simultaneously provide power to, and 
measure the resistance response of the heater element, the heater element is preferably provided in one leg of a 
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Wheatstone bridge, for example as shown in Figure 10. Thus, in the illustrative embodiment, the heater block 510 
preferably includes a circuit similar to that shown in Figure 10. 

[0061] The modulated heater input signal is provided to a power input terminal of the Wheatstone bridge, such as 
power input terminal 354 of Figure 10, which is either directly or indirectly coupled to the heater element. In this con- 
5 figuration, the Wheatstone bridge provides a differential output signal with an amplitude that is proportional to the 
resistance of the heater element. The differential output of the Wheatstone bridge may be provided to a differential 
amplifier, as shown in Figure 1 0, or may be directly provided as the output of heater block 51 0. An illustrative transient 
resistance response of the heater element is shown in Figure 15. 

[0062] For accurate measurement of the heater resistance, the high frequency component may be removed from 
10 the heater output signal of the Wheatstone bridge using any number of techniques. One such technique is to provide 
a low-pass filter 512 at the output of the heater block 510 which allows only the lower frequency component of the 
heater output signal to pass. The resulting filtered signal may then be provided to a stop input of a first high frequency 
timer 514. Preferably, the heater output signal of the heater block 510 is AC coupled to the filter 512, as shown by 
capacitor 516. The filter 512 may provide the AC coupling function, or a separate element such as capacitor 516 may 
15 be provided. 

[0063] The start input of the first high frequency timer 514 may be coupled to the lower frequency component of the 
heater input signal. More preferably however, the lower frequency component of the heater input signal is coupled to 
the start input of the high frequency timer 514 via a bulk resistor 520 and a capacitor 522 as shown. In any case the 
start input of the high frequency timer 514 is responsive to the lower frequency component of the heater input signal. 

20 The bulk resistor 520 preferably has a low temperature coefficient to minimize any phase lag caused thereby. Capacitor 
522 AC couples the lower frequency component of the input signal to the start input of the high frequency timer 514. 
[0064] Because both the heater input signal and the heater output signal are preferably AC coupled to the timer start 
and stop inputs, respectively, of the timer 514, the zero crossing points of the heater input and output signals may be 
used to trigger the high frequency timer 514. This allows the time lag measurement of the high frequency timer 514 to 

25 be relatively independent of the amplitudes of the input and output signals, thereby increasing the accuracy of the 
measurement. 

[0065] As can be seen from the above description, high frequency timer 514 starts when the AC coupled lower 
frequency component of the heater input signal crosses some predefined threshold, preferably zero. Likewise, high 
frequency timer 514 stops when the AC coupled heater output signal, which represents the resistance of the heater 
element, crosses some predefined threshold, preferably zero. The result is a heater time lag Az h . When the heater is 
exposed to a fluid at substantially zero flow, the heater time lag Az h equals Az k p . As described below, the thermal 
conductivity of the fluid can be determined from Az k p . 

[0066] In some applications, it may be desirable to control the amplitude of the transient elevated temperature con- 
dition of the heater element. This is preferably accomplished by providing an amplitude control signal that is indicative, 
of the amplitude of the resistance change in the heater element. The amplitude control signal may be provided by a 
rectifier 526, which rectifies the filtered output signal as shown. The generator 500 may include an amplitude control 
block 530, which accept the amplitude control signal via interface 528, and adjusts the amplitude of the lower frequency 
component and the inverted copy thereof such that the amplitude of the resistance change in the heater element 
remains at a relatively constant level. 
*o [0067] With respect to the first sensor element, a second generator 550 generates and provides the high frequency 
component of the sensor input signal to a modulator 552 via interface 554. The generator 550 also generates the lower 
frequency component, and an inverted copy thereof, and provides these signals to modulator 552 via interfaces 556 
and 558, respectively. As described above, the modulator 552 modulates the high frequency component using the 
lower frequency component signals to produce a modulated sensor input signal as shown in Figure 13. The power 
4 $ delivered by the modulated sensor input signal is shown in Figure 14. 

[0068] The modulated sensor input signal is provided to sensor block 560. To simultaneously provide power to : and 
measure the resistance response of the sensor element, the second sensor element is preferably provided in one leg 
of a Wheatstone bridge : for example as shown in Figure 10. Thus, in the illustrative embodiment, the sensor block 560 
preferably includes a circuit similar to that shown in Figure 10. 
so [0069] The modulated sensor input signal is provided to a power input terminal of the Wheatstone bridge, such as 
power input terminal 380 of Figure 1 0, which is either directly or indirectly coupled to the sensor element. In this con- 
figuration, the Wheatstone bridge provides a differential output signal with an amplitude that is proportional to the 
resistance of the sensor element. The differential output of the Wheatstone bridge may be provided to a differential 
amplifier, as shown in Figure 1 0, or may be directly provided as the output of sensor block 560. An illustrative transient 
55 resistance response of the sensor element is shown in Figure 15. 

[0070] For accurate measurement of the sensor resistance, the high frequency component of the modulated heater 
input signal may be removed using any number of techniques. One such technique is to provide a low-pass filter 562 
at the output of the sensor block 560 which only allows the lower frequency component of the sensor output signal to 
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pass. The resulting filtered signal may then be provided to a stop input of a second high frequency timer 564. Preferably, 
the sensor output signal of the sensor block 560 is AC coupled to the filter 562, as shown by capacitor 566. The filter 
562 may provide the AC coupling function, or a separate element such as capacitor 566 may be provided. 
[0071] The start input of high frequency timer 564 may be coupled to the lower frequency component of the sensor 

5 input signal. More preferably, however, the lower frequency component of the heater input signal is coupled to the start 
input of the high frequency timer 564 via a bulk resistor 570 and a capacitor 572 as shown. The bulk resistor 570 
preferably has a low temperature coefficient to minimize any phase lag caused thereby. Capacitor 572 AC couples the 
lower frequency component of the sensor input signal to the start input of the high frequency timer 564. 
[0072] Because both the sensor input signal and the sensor output signal are AC coupled to the timer start and stop 

10 inputs, respectively, of the timer 564, the zero crossing points of the sensor input and output signals may be used to 
trigger the high frequency timer 564. This allows the time lag measurement of the high frequency timer 564 to be 
relatively independent of the amplitudes of the sensor input and output signals, thereby increasing the accuracy of the 
measurement. 

[0073] As can be seen from the above description, high frequency timer 564 starts when the AC coupled lower 
is frequency component of the sensor input signal crosses some predefined threshold, preferably zero. Likewise, high 
frequency timer 564 stops when the AC coupled sensor output signal, which represents the resistance of the second 
sensor element, crosses some predefined threshold, preferably zero. The result is a first sensor time lag Az s1 . 
[0074] To provide a sensor input signal that has the proper frequency, phase and amplitude to cause the resulting 
temperature response in the first sensor element to track the temperature disturbance of the fluid, a phase shifter 580 
20 and analog divider 582 may be provided. The desired phase and amplitude of the sensor input signal can be determined 
during a calibration procedure. In a preferred embodiment, the sensor element is first subjected to a vacuum condition, 
and a sensor input signal is provided to the first sensor element as described above. Because no fluid surrounds the 
sensor element, substantially no heat is transferred from the sensor element to the fluid. A sensor time lag Az s1vac is 
measured by the high frequency timer 564, and subsequently stored. An illustrative sensor time lag Az s1vac is shown 
25 jn Figure 16, and represents the delay between the sensor input signal 551 and the resistive response 553 of the 
sensor element. In addition, a power/resistance ratio between the sensor input signal and the resulting resistance 
response of the sensor element is measured and stored. 

[0075] The sensor element is then subjected to a fluid of interest. During operation, the phase of the sensor input 
signal is adjusted by phase shifter 580 so that the resulting first sensor time lag Az s1> measured by high frequency 

30 timer 564, equals the sensor time lag measured under vacuum conditions Az s1vac . Generator 550 may include a com- 
parator, processor or other means for comparing the measured time lag Az s1 and the calibrated time lag Az s1vac . 
[0076] Likewise, the power of the input signal (P) and the resulting resistance response (R) of the first sensor element 
may be provided to analog divider 582, as shown. Analog divider 582 may divide these signals to provide a power/ 
resistance ratio. Generator 550 may then adjust the amplitude of the first sensor input signal so that the resulting power/ 

35 resistance ratio equals the power/resistance ratio recorded under vacuum conditions. 

[0077] Accordingly, the first sensor input signal heats the first sensor element in phase and with the proper amplitude 
to match the temperature disturbance in the fluid such that substantially no heat is transferred from the fluid to the first 
sensor element. This may help reduce or eliminate the first sensor time lag Az s1 as a factor when determining the first 
transit time Az 1 of the temperature disturbance. 

^0 [0078] Timer 590 directly measures the transit time Az 1 of the temperature disturbance from the heater element to 
the first sensor element. Timer 590 is started when the response of the heater element crosses a predefined threshold, 
preferably zero. Timer 590 is stopped with the response of the first sensor element crosses a predefined threshold, 
preferably zero. Since the first sensor is tracking the temperature disturbance in the fluid, the first sensor time lag Az s1 
is removed as an error source. When the first sensor is exposed to the fluid at substantially zero flow, the first sensor 

45 time lag Az s1 can be used to determine the thermal diffusivity and specific heat of the fluid, as further described below. 
[0079] Finally, and with respect to the second sensor element, a third generator 600 generates and provides the high 
frequency component of the second sensor input signal to a modulator 602 via interface 604. The generator 600 also 
generates the lower frequency component, and an inverted copy thereof, and provides these signals to modulator 602 
via interfaces 606 and 608, respectively. As described above, the modulator 602 modulates the high frequency com- 

50 ponent using the lower frequency component signals to produce a modulated sensor input signal as shown in Figure 
13. The power delivered by the modulated sensor input signal is shown in Figure 14. 

[0080] The modulated sensor input signal is provided to sensor block 61 0. To simultaneously provide power to : and 
measure the resistance response of the sensor element, the sensor element is preferably provided in one leg of a 
Wheatstone bridge, for example as shown in Figure 10. Thus, in the illustrative embodiment, the sensor block 610 
55 preferably includes a circuit similar to that shown in Figure 10. 

[0081] The modulated sensor input signal is provided to a power input terminal of the Wheatstone bridge, such as 
power input terminal 401 of Figure 10, which is either directly or indirectly coupled to the sensor element. In this con- 
figuration, the Wheatstone bridge provides a differential output signal with an amplitude that is proportional to the 
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resistance of the sensor element. The differential output of the Wheatstone bridge may be provided to a differential 
amplifier, as shown in Figure 1 0, or may be directly provided as the output of sensor block 61 0. An illustrative transient 
resistance response of the sensor element is shown in Figure 15. 

[0082] For accurate measurement of the sensor resistance, the high frequency component of the modulated heater 
5 input signal may be removed using any number of techniques. One such technique is to provide a low-pass filter 612 
at the output of the sensor block 610 which only allows the lower frequency component of the sensor output signal to 
pass. The resulting filtered signal may then be provided to a stop input of a second high frequency timer 61 4. Preferably, 
the sensor output signal of the sensor block 610 is AC coupled to the filter 612, as shown by capacitor 616. The filter 
612 may provide the AC coupling function, or a separate element such as capacitor 61 6 may be provided. 
10 [0083] The start input of the high frequency timer 614 may be coupled to the lower frequency component of the 
sensor input signal. More preferably, however, the lower frequency component of the sensor input signal is coupled to 
the start input of the high frequency timer 614 via a bulk resistor 620 and a capacitor 622 as shown. The bulk resistor 
620 preferably has a low temperature coefficient to minimize any phase lag caused thereby. Capacitor 622 AC couples 
the lower frequency component of the sensor input signal to the start input of the high frequency timer 614. 
[0084] Because both the sensor input signal and the sensor output signal are AC coupled to the timer start and stop 
inputs, respectively, of the timer 614, the zero crossing points of the sensor input and output signals may be used to 
trigger the high frequency timer 614. This allows the time lag measurement of the high frequency timer 614 to be 
relatively independent of the amplitudes of the sensor input and output signals, thereby increasing the accuracy of the 
measurement. 

20 [0085] As can be seen from the above description, high frequency timer 614 starts when the AC coupled lower 
frequency component of the third sensor input signal crosses some predefined threshold, preferably zero. Likewise, 
high frequency timer 61 4 stops when the AC coupled sensor output signal, which represents the resistance of the third 
sensor element, crosses some predefined threshold, preferably zero. The result is a second sensor time lag Az^. 
[0086] To provide a sensor input signal that has the proper frequency, phase and amplitude to cause the resulting 

25 temperature response in the second sensor element to track the temperature disturbance of the fluid, a phase shifter 
630 and analog divider 632 may be provided. The desired phase and amplitude of the sensor input signal can be 
determined during a calibration procedure. In a preferred embodiment, the second sensor element is first subjected 
to a vacuum condition, and a sensor input signal is provided to the second sensor element as described above. Because 
no fluid surrounds the sensor element, substantially no heat is transferred from the sensor element to the fluid. A 

30 second sensor time lag Az s2vac is measured by the high frequency timer 614, and subsequently stored. An Illustrative 
second sensor time lag Az^^ is shown in Figure 16, and represents the delay between the sensor input signal 551 
and the resistive response 553 of the second sensor element. In addition, a power/resistance ratio between the sensor 
input signal and the resulting resistance response of the second sensor element is measured and stored. 
[0087] The third sensor element is then subjected to a fluid of interest. During operation, the phase of the sensor 

35 input signal is adjusted by phase shifter 630 so that the resulting second sensor time lag Az s2 , measured by high 
frequency timer 61 4 ; equals the second sensor time lag measured under vacuum conditions Az s2vac . Generator 600- 
may include a comparator, processor or other means for comparing the measured time lag Az s2 and the calibrated 
time lag Az s2vac . 

[0088] Likewise, the power of the input signal (P) and the resulting resistance response (R) of the second sensor 
40 element may be provided to analog divider 632, as shown. Analog divider 632 may divide these signals to provide a 
power/resistance ratio. Generator 600 may then adjust the amplitude of the second sensor input signal so that the 
resulting power/resistance ratio equals the power/resistance ratio recorded under vacuum conditions. 
[0089] Accordingly, the second sensor input signal heats the second sensor element in phase and with the proper 
amplitude to match the temperature disturbance in the fluid such that substantially no heat is transferred from the fluid 
to the second sensor element. This may help reduce or eliminate the second sensor time lag Az^ as a potential error 
source when determining the transit time of the temperature disturbance. 

[0090] Timer 650 direclly measures Ihe transit lime Az 2 of the temperature disturbance from the heater element to 
the second sensor element. Timer 650 is started when the response of the heater element crosses a predefined thresh- 
old, preferably zero. Timer 650 is stopped with the response of the second sensor element crosses a predefined thresh- 
old, preferably zero. Since the second sensor is tracking the temperature disturbance in the fluid, the sensor time lag 
Az s2 is eliminated as an error source. 

[0091] The velocity of the fluid can be determined from the transit time values Az 1 and Az 2 using the abovc-dcscribcd 
relation: 



45 



50 



55 



v 2 = (d 1 2 /z 1 -d 2 2 /z 2 )/( Zl -z 2 ) 



(8) 
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[0092] Thus, the illustrative embodiment of Figure 11 A and Figure 11 B can be used to calculate the velocity of the 
fluid of interest relatively independently of the properties of the fluid, while minimizing the effects of the heater and 
sensor time lags. 

[0093] Further, the illustrative embodiment of Figure 11 A and Figure 11 B can be used to determine the thermal 
conductivity, thermal diffusivity and specific heat of the fluid of interest, as described more fully below. 

THERMAL CONDUCTIVITY 

[0094] The thermal dynamic interaction between a heater element that is forced to experience fluctuating tempera- 
tures of: 



T = T 0 + T.,sin(u>z +T) 



(9) 



15 in response to a fluctuating input power 



20 



p = P Q (1 + sin(coz)) 
can be described by a simple differential equation 



(10) 



c pv t dT/dz 



P 0 (1 + sin(u>z)) - (h, + h 2 )(T- T f ) - h 3 (T-T b ) 



(11) 



25 



30 



35 



40 



45 



50 



55 



TABLE I 



Symbol 



Nomenclature 



f 

t 

T 



h 2 
h 3 
Li 
L 2 

z 
T 

AZ, 

Az h 
AZ S 

AZ 1 
AZ 2 

di 
d 2 



k.p 



frequency of the input signal, H z 
2jcf, H 2 

specific heat per unit volume for the heater film and support member (1 0% Platinum, 90% S13N4 
Microbridge composite, j/(cm 3 k) 
heater film thickness, cm 

sensor base temperature, with peak-to-peak amplitude of 2T 0 , k 
fluid temperature, k 
substrate temperature, k 

coefficient of conductive heat transfer to the fluid of interest = k/L,), W/cm 3 

coefficient of forced convective heat transfer to the fluid interest under laminar flow (= k/l^), W/cm 3 
coefficient of conductive heat transfer to the substrate, W/cm 3 

characteristic length of thermal conduction from the heater element into the fluid phase, cm 
characteristic length of convective heat transfer, cm 
time, s 
Temperature 

time lag between the heater input signal and the resistance response of the heater means when 
measured with substantially zero fluid flow 

time lag between the heater input signal and the resistance response of the heater means when 
measured with fluid flow 

time lag between the sensor input signal and the resistance response of the corresponding sensor means 
when measured with fluid flow 

time lag between the sensor input signal and the resistance response of the sensor means under vacuum 
conditions 

time lag required for a thermal disturbance to travel between the heater and the first sensor 
time lag required for a thermal disturbance to travel between the heater and the second sensor 
separation distance between the heater element and the first sensor element 
separation distance between the heater element and the second sensor element (d1 1 62) 
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15 



20 



25 



30 



35 



45 



55 



Symbol 



Nomenclature 



phase lag between input signal and the resistance of the heater means (T= Az27if), radians 



Integration of equation (11) leads to the solution for the phase lag, Y. and the DC and AC signal amplitudes, T 0 and 
T 1t respectively as follows: 

io T = arctan(-27ifc pv t/(h 1 + h 2 + h 3 )) ' (12) 

Az = T/(27if) (13) 
T 0 = ((h 1 + h 2 )T f + h 3 T b + P 0 )/(h 1 + h 2 + h 3 ) (14) 



T 1 - P 0 /((h 1 + h 2 + h 3 ) 2 + (c^to)) 2 )* (15) 

[0095] The contributions of h 1} h 2 and h 3 to the phase lag T can be isolated and individually measured. During a 
calibration procedure, for example, the value of h 3 can be determined by subjecting the heater element to a vacuum 
condition, thereby reducing h t and h 2 to zero. A time lag value may then be measured between the input signal and 
the output signal of the heater element under the vacuum condition. The value of h 2 may then be calculated using the 
relation: 

h 3 = -27ifc pv t/tan(T) (16) 

[0096] The value of may then be determined by subjecting the heater element to the fluid of interest at atmospheric 
pressure and substantially zero flow, thereby reducing h 2 to zero. A time lag can then be measured between the input 
signal and the output signal of the heater element under atmospheric pressure. The value of h 1 can then be calculated 
using the relation: 



h 1 = [-2jifc pv tAan(T)] - h 3 (17) 

where h 3 is known from above. 
40 [0097] Finally, the value of h 2 may be determined by subjecting the heater means to the fluid of interest at a prede- 
termined non-zero flow rate. A time lag can then be measured between the input signal and the output signal of the 
heater means under the non-zero flow condition. The value of h 2 can then be calculated using the relation: 



h 2 = l-27cfc pv tAan(T)j - h, - h 3 (18) 



where h 2 and h 3 are known from above. 

[0098] In the illustrative embodiment shown in, for example, Figure 5, the heater element 224 and the support member 
230 have a corresponding specific heat value, C pv . Further, the heater element 224 has a coefficient of conductive 
50 heat transfer to the substrate 232, h 3 . Once these parameters are determined, for example by prior calibration as 
described above, the thermal conductivity, k, of the fluid of interest can be determined using the relation: 

k = (-27cfc pv t/tan(Az k p 2nf) - h 3 )L 1 (19) 
THERMAL PIFFUSIVITY AND SPECIFIC HEAT 

[0099] The propagation or spread of a local, non-uniform temperature distribution can be derived from the following 
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three-dimensional temperature diffusivity or heat conductivity equation, to which a term has been added to represent 
fluid velocity in the x-direction: 

5 aT/dz = D x d 2 T/dx 2 - v x 8T/dx (20) 

The one dimensional transient solution to equation (20) is: 

w T ~~ z* a5 exp(-(d 1 - vz) 2 /(4D t z)) (21) 

where x has been replaced by (± d 1 ± vz) to include the velocity term. To determine the time lag between very short 
forcing functions at time z=0, and the arrival of these pulses in terms of maxima in T(z) at the first sensor element 
position d 1t equation (21) is differentiated and set to zero: 



15 



20 



25 



cff/dz = 0 = -1/z 2 i (1/z)(2v(d 1 - vz)/4D t z + (d 1 - vz) 2 /4D t z 2 ) (22) 



multiplying by z 2 and rearranging yields the relation: 



-4D l z-v 2 z 2 + d 1 2 = 0 (23) 

Because the thermal conductivity, thermal diffusivity and specific heat of the fluid of interest are preferably determined 
at substantially zero flow, v is set to zero, resulting in the following relation for D t : 

D^d^Az,) (24) 

30 where z has been replaced by z 1 . which represents the transit time for the temperature disturbance to travel a distance 
d 1) or in this case, from the heater element to the first sensor element. 
[0100] The specific heat of the fluid of interest may then be determined using the relation: 



35 c p = k/D t (25) 

where the thermal conductivity k, and the thermal diffusivity, D t , are determined as described above. 

[0101] Figure 17 is a timing diagram showing the various time lags measured by the embodiment of Figure 11 A 

and Figure 11B. The heater time lag value measured with substantially zero flow Az kp is determined as described 

above, and represents the time lag between the heater input signal shown at 700 and the response of the heater 

element shown at 702. The heater time lag Az k p is preferably used to calculate the thermal conductivity of the fluid of 

interest. 

[01 02] The sensor time lags Az s1 and Az s2 are preferably forced to match the sensor time lag values measured under 
vacuum conditions Az s1 2vac . Thus, the sensor input signal shown at 708 is provided with the proper phase relative to 
the heater input signal, and with the proper amplitude, so that the substantially no heat is transferred from the fluid to 
the corresponding sensor element. The temperature response of the corresponding sensor element is shown at 708. 
[0103] The first transit time AZ-, is the time lag between the heater response and the sensor response of the first 
sensor, preferably with Az s1 = Az s1vac as shown. As described above, the transit time Az^ can be used to calculate the 
thermal diffusivity and the specific heat of the fluid of interest, when measured with substantially zero fluid flow. The 
transit time AZ 2 is the time lag between the heater response and the sensor response of the second sensor, preferably 
with Az s2 = Az s2vac as shown. The transit time Az 2 can also be used to calculate the thermal diffusivity and the specific 
heat of the fluid of interest when measured with substantially zero fluid flow. Finally, the transit times Az 1 and Az 2 can 
be used to determine the velocity of the fluid with f!ow ; relatively independently of the properties of the fluid. 
[0104] It is contemplated that the heater and sensor input signals may be active during a measurement period, and 
inactive during a wail period. Further, it is contemplated that the wail period may be longer than the measurement 
period. This may conserve power, which may be particularly important when the sensor system is battery powered. 
[0105] Figure 18 is a block diagram showing a processing means that is used for calculating the velocity, thermal 
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conductivity and specific heat of the fluid of interest using the time lag values provided by Figure 11 A and Figure 1 1 B. 
Using the above-described relations, processing means 800 may accept Az kp , f, c pv , t, h3 and L 1 to calculate the 
thermal conductivity, k, of the fluid of interest. Likewise, processing means 800 may use the thermal conductivity, k, 
one of the transit times such as Az 1 when measured with substantially zero fluid flow, and the distance between the 
heater and the corresponding sensor "d M to determine the thermal diffusivity and specific heat of the fluid of interest. 
[0106] Finally, processing means 800 may use the first and second transit times Az 1 and Az 2 when measured with 
fluid flow, along with the distances between the heater element and the corresponding sensor elements d 1 and d 2 , to 
determine the velocity of the fluid. 

[0107] Figure 19 and Figure 20 are schematic diagrams of illustrative embodiments that use optical sensor means 
to measure the transit time required for a thermal disturbance to travel from a first location to a second location, In both 
embodiments, a heater input signal is provided to a heater element 900 to induce thermal disturbances in the fluid. In 
Figure 19, a first light source 902 and a second light source 904 are provided at a first and a second distance from 
the heater element 900, respectively. A corresponding first light sensor 906 and second light sensor 908 are provided, 
whereby each light sensor is separated from the corresponding light source so that the fluid passes therebetween. 
Each of the light source/sensor pairs detects the thermal disturbance by measuring a deflection of the light source due 
to the modulated index of refraction of the heated fluid. This provides a temperature sensor with little or no thermal 
mass, thereby minimizing the effects of the sensor time lags inherent in temperature sensitive resistive type sensors. 
[0108] In Figure 20, two light source/sensor pairs 916 and 918 are provided adjacent to one another. The light is 
emitted by each or the light sources, reflected from a reflective surface 920, and returned to the corresponding light 
sensor. As in Figure 1 9, each of the light source/sensor pairs detects the thermal disturbance by measuring a deflection 
of the light source due to the modulated index of refraction of the heated fluid. 

[0109] In both embodiments, the light source may any type of light source, but preferably is a LED (Light Emitting 
Diode) or Vertical Cavity Surface Emitting Laser (VCSEL). The sensor may be any type of light sensor, but preferably 
a photodiodc (PD) or RCPD (Resonant Cavity PhotoDiode). 

[0110] Having thus described the preferred embodiments of the present invention, those of skill in the art will readily 
appreciate that the teachings found herein may be applied to yet other embodiments within the scope of the claims 
hereto attached. 



Claims 

1. Apparatus for determining the fluid velocity, v, of a fluid of interest, comprising: 

heater means in thermal communication with the fluid of interest, said heater means having a resistance that 
changes with temperature: 

heater energizing means connected to said heater means for energizing said heater means, said heater en-, 
ergizing means providing a periodic time-varying input signal to said heater means to induce a transient ele- 
vated temperature condition in said heater means; 

first output means for providing a first output signal that is proportional to the resistance of said heater means; 
sensor means in thermal communication with the fluid of interest, the sensor means having a resistance that 
changes with temperature; 

sensor energizing means for energizing the sensor means : said sensor energizing means providing a periodic 
time-varying input signal to the sensor means that is out of phase relative to the periodic time-varying input 
signal of said heater energizing means, whereby the sensor input signal has the same frequency as the heater 
input signal; 

sensor output means for providing a sensor output signal that is proportional to the resistance of the sensor 

means; 

first time lag means for determining a first time lag value corresponding to the delay between generating a 
transient elevated temperature condition at said heater means and sensing a transient elevated temperature 
response at said sensor means; 

second time lag means for determining a second time lag value representing the delay between corresponding 
points of the sensor input signal provided by the sensor energizing means and the sensor output signal, 
phase shifter means for shifting the phase between the periodic time-varying input signal provided by said 
heater energizing means and the periodic time-varying input signal provided by the sensor energizing means 
so that the second time lag value substantially equals a calibration second time lag value measured when the 
sensor means is exposed to a vacuum; and 

determining means for determining the fluid velocity, v, of the fluid of interest using the first time lag value. 
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2. Apparatus according to claim 1 wherein the sensor input signal provided by said sensor energizing means has an 
amplitude that is selected to provide a power-resistance ratio in the sensor means that is substantially identical to 
the power-resistance ratio in the sensor means when the sensor is exposed to a vacuum. 

5 3. A method for determining the fluid velocity, v, of a fluid of interest using a heater element and at least two spaced 
sensor elements, the method comprising the steps of: 

providing a first one of the sensor elements a first distance from the heater element, and a second one of the 
sensor elements a second distance from the heater element, wherein the second distance is greater than the 
to first distance; 

energizing said heater element to providing a transient elevated temperature condition in said heater element; 
determining a first time lag value corresponding to the time delay between generating the heater input signal 
and sensing the first sensor output signal; 

determining a second time lag value between generating the heater input signal and sensing the second sensor 
is output signal; 

determining the fluid velocity, v, of the fluid of interest using the first and second time lag values and using the 
relation: 

20 v = {(d 1 2 /Az 1 - d 2 2 /Az 2 )/( A Zl -Az 2 )}° 5 

where, 

d 1 = the distance between the heater element and a first sensor element; 

25 

d 2 - the distance between the heater element and the second sensor element, where Id^ ld 2 l 

Az t = the first time lag value; and 

^0 Az 2 = the second time lag value. 

4. A method according to claim 3 further comprising the steps of: 

energizing said first sensor element with a first sensor input signal having a phase and amplitude such that 
35 substantially no heat is transferred from the fluid of interest to the first sensor element during the transient 

elevated temperature condition; and 

energizing said second sensor element with a second sensor input signal having a phase and amplitude such 
that substantially no heat is transferred from the fluid of interest to the second sensor element during the 
transient elevated temperature condition. 



40 



5. Apparatus for determining the fluid velocity, v, of a fluid of interest, comprising: 



heater means in thermal communication with the fluid of interest, said heater means having a resistance that 
changes with temperature; 

45 heater energizing means connected to said heater means for energizing said heater means, said heater en- 

ergizing means providing a transient elevated temperature condition in said heater means; 
at least two sensor means in thermal communication with the fluid of interest, each of the at least two sensor 
means spaced a different distance from said heater means, and each having a resistance that changes with 
temperature; 

50 time lag means for determining at least two time lag values each corresponding to the delay between said 

heater means and a corresponding sensor means; and 

determining means for determining the fluid velocity, v, of the fluid of interest using the at least two time lag 
values. 

55 6. A method for determining the fluid velocity, v, of a fluid of interest using a heater element and at least two spaced 
sensor elements, the method comprising the steps of: 

providing a first one of the sensor elements a first distance from the heater element, and a second one of the 
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sensor elements a second distance from the heater element, wherein the second distance is greater than the 
first distance; 

energizing said heater element to providing a transient elevated temperature condition in said heater element: 
determining a first time lag value between the heater element and the first one of the sensor elements; 
determining a second time lag value between the heater element and the second one of said sensor elements; 
determining the fluid velocity, v, of the fluid of interest using the first and second time lag values. 
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Fig. 1 (PRIOR ART) 




Fig. 2 (PRIOR ART) 




Fig. 3 (prior art) 
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Fig. 6 
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